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The  computational  modeling  and  design  of  an  actively-cooled  microvascular  fin  specimen  is  presented. 
The  design  study  is  based  on  three  objective  functions:  (i)  minimizing  the  maximum  temperature  in 
the  thermally  loaded  fin,  (ii)  optimizing  the  flow  efficiency  of  the  embedded  microchannel,  and  (iii)  min¬ 
imizing  the  void  volume  fraction  of  the  microvascular  material.  A  recently  introduced  Interface-enriched 
Generalized  Finite  Element  Method  (IGFEM)  is  employed  to  evaluate  the  temperature  field  in  a  2D  model 
of  the  specimen,  allowing  for  the  accurate  and  efficient  capturing  of  the  gradient  discontinuity  along  the 
fluid/solid  interface  without  the  need  of  meshes  that  conform  to  the  geometry  of  the  problem.  Finding  the 
optimal  shape  of  the  embedded  microchannel  is  thus  accomplished  with  a  single  non-conforming  mesh 
for  all  configurations.  Prior  to  the  optimization  study,  the  IGFEM  solver  is  validated  through  comparison 
with  infrared  measurements  of  the  thermal  response  of  an  epoxy  fin  with  a  sinusoidal  microchannel. 
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1.  Introduction 

Modeled  after  circulatory  systems  found  in  many  living  organ¬ 
isms,  microvascular  materials  with  microchannels  of  diameters 
ranging  from  lOjim  to  1  mm,  embedded  in  polymeric  matrices, 
have  shown  great  promise  in  a  variety  of  engineering  applications 
such  as  autonomic  healing  and  active  cooling  [1-3].  Recent 
advances  in  manufacturing  techniques  for  microfluidic  materials 
and  devices  [4,5]  have  extended  the  application  of  such  systems 
to  a  broad  field  including  autonomic  materials  [6],  biotechnology 
[7,8],  chemical  reactors  [9],  and  microelectromechanical  systems 
(MEMS)  [10-12].  In  many  of  these  applications,  biomimicry  has 
been  used  as  an  inspiration  for  the  design  of  the  microvascular  sys¬ 
tem,  while  manufacturing  constraints  often  play  a  key  role  in 
restricting  the  design  space.  Nevertheless,  the  high  number  of 
parameters  defining  the  optimal  configuration  of  the  embedded 
microchannels  topology,  thermal  loading,  etc.,  make  the  numerical 
design  of  such  materials  quite  challenging. 

The  focus  of  the  present  study  is  on  the  computational  design  of 
an  actively-cooled  microvascular  fin.  The  circulation  of  a  coolant  in 
microchannels  reduces  the  temperature  of  the  surrounding 
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material  in  two  specific  ways:  (i)  by  direct  extraction  of  heat  from 
the  thermally  loaded  material,  and  (ii)  by  redistributing  the  heat 
between  hotter  and  colder  regions  of  the  domain.  The  first  mecha¬ 
nism  is  effective  if  the  length  of  the  embedded  microchannels 
between  the  inlet  and  outlet  is  rather  short,  i.e.,  when  the  fluid  ex¬ 
its  before  it  becomes  too  hot  and  looses  its  cooling  capacity.  MEMS 
with  an  active  cooling  substrate  is  an  example  of  this  type  of 
microvascular  materials  [12,10].  For  microvascular  materials  with 
larger  dimensions  and  consequently  larger  distances  between  the 
inlet  and  outlet  of  the  channels,  removing  the  heat  from  the 
domain  through  active  cooling  may  not  effectively  reduce  the  tem¬ 
perature  of  the  domain.  In  this  case,  the  extracted  thermal  energy 
through  the  fluid  only  accounts  for  a  small  portion  of  the  thermal 
energy  given  to  the  system.  Instead,  the  microvascular  network 
must  be  designed  to  effectively  redistribute  the  heat  inside  the 
material  by  extracting  the  heat  from  warmer  areas  of  the  domain 
and  releasing  it  in  colder  regions.  An  example  of  nature’s  designs 
for  this  type  of  actively-cooled  systems  is  the  arrangement  of  blood 
vessels  for  the  thermal  regulation  in  human  skin  [13,14].  Depend¬ 
ing  on  environmental  conditions,  the  vascular  response  of  the  skin 
is  able  to  adjust  its  temperature  to  preserve  or  dissipate  the  heat 
when  the  outside  temperature  is  decreased  or  increased,  respec¬ 
tively  (Fig.  1(a)). 

Additional  motivation  for  the  work  presented  hereafter  stems 
from  recent  advances  in  the  manufacturing  of  large,  microvascular 
3D  woven  composites  based  on  a  novel  approach  with  sacrificial 
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Fig.  1.  Biological  and  biomimetic  examples  of  microvascular  systems  for  active 
cooling:  (a)  microvascular  blood  vessels  for  thermal  regulation  in  the  human  skin 
(inspired  from  [13]);  (b)  micro-CT  image  of  sinusoidal-shaped  microchannels  (filled 
with  eutectic  Ga-In  alloy  for  a  higher  resolution)  embedded  in  a  3D  woven  glass 
fiber/epoxy  matrix  composite  created  with  the  sacrificial  fiber  manufacturing 
technique. 


fibers  [5].  This  method  allows  for  the  creation  of  3D  glass  fiber/ 
epoxy  matrix  composites  with  sinusoidal  microchannels  fully  inte¬ 
grated  into  the  material  microstructure  (Fig.  1(b)). 

Inspired  by  these  two  examples,  we  adopt  a  sinusoidal  shape  for 
the  embedded  microchannel  in  the  design  study  of  a  microvascular 
actively-cooled  fin,  and  aim  to  determine  the  optimal  values  of  the 
amplitude  and  wavelength  of  the  microchannel.  While  the  main 
objective  is  to  evaluate  the  optimal  geometry  of  the  microchannels 
to  minimize  the  maximum  temperature  of  the  domain,  their 
impact  on  mechanical  properties  such  as  the  stiffness  of  the  micro¬ 
vascular  material  is  also  addressed  indirectly  by  minimizing  the 
void  volume  fraction  associated  with  the  embedded  network.  The 
third  objective  function  adopted  in  the  design  study  is  the  minimi¬ 
zation  of  the  pressure  drop  needed  to  circulate  the  coolant,  i.e.,  the 
flow  efficiency  of  the  microchannels.  Similar  optimization  prob¬ 
lems  for  the  design  of  microvascular  materials  have  been 
addressed  in  [2,15-19]. 

Although  the  geometrical  setting  of  the  2D  shape  optimization 
study  is  relatively  simple,  the  efficiency  and  accuracy  of  the  finite 
element  solver  used  to  predict  the  impact  of  the  sinusoidal  micro- 
channel  on  the  thermal  response  of  the  microvascular  fin  is  a  crit¬ 
ical  aspect  of  this  study.  Employing  the  standard  finite  element 
method  (FEM)  for  such  problems  is  quite  cumbersome  since  this 
method  requires  for  each  configuration  the  creation  of  a  mesh  that 
conforms  to  the  microstructure  of  the  material  to  capture  the  dis¬ 
continuity  in  the  gradient  field  along  the  fluid/solid  interface.  In¬ 
stead,  we  adopt  in  this  paper  a  recently  introduced  generalized 
finite  element  method  [20,21]  to  evaluate  the  temperature  field 
in  the  microvascular  fin.  As  in  conventional  Generalized  Finite  Ele¬ 
ment  Methods  (GFEM)  [22-24],  an  a  priori  knowledge  of  the  solu¬ 
tion  field  in  the  form  of  weakly  discontinuous  enrichment 
functions  is  incorporated  in  the  numerical  approximation  of  the 
thermal  field.  However,  instead  of  introducing  nodal  enrichment 
degrees  of  freedom  as  in  conventional  GFEM,  the  interface-based 
GFEM  formulation  (referred  to  hereafter  as  IGFEM)  applies  enrich¬ 
ments  at  the  nodes  created  by  the  intersection  of  the  interface  with 
edges  of  the  non-conforming  element. 


The  outline  of  the  paper  is  as  follows:  In  Section  2,  we  first 
present  the  strong  and  weak  forms  of  the  governing  equations 
describing  the  temperature  field  in  the  actively  cooled  microvascu¬ 
lar  material  of  interest.  We  then  introduce  the  IGFEM  formulation 
for  this  problem  and  briefly  discuss  some  of  the  key  characteristics 
of  the  numerical  scheme.  In  Section  3,  we  summarize  the  results 
of  a  validation  study  for  the  IGFEM  solver,  based  on  direct  compar¬ 
isons  with  experimental  observations  of  the  thermal  field  in  a 
microvascular  epoxy  fin.  Finally,  we  perform  in  Section  4  a  shape 
optimization  study  for  a  sinusoidal  microchannel  embedded  in  a 
long  polymeric  fin  subjected  to  a  uniform  heat  flux  along  one  of  its 
faces,  and  investigate  how  the  predicted  optimal  shape  of  the  micro- 
channel  is  affected  by  the  choice  of  thermal  boundary  conditions. 

2.  Governing  equations  and  IGFEM  formulation 

To  model  the  temperature  field  in  an  actively-cooled  microvas¬ 
cular  material,  we  solve  the  convection-diffusion  equations  on  a 
2D  domain  Q  C  R2  composed  of  two  mutually  exclusive  regions 
Qf  and  Qs ,  corresponding  to  the  fluid  and  solid  phases  in  a  micro¬ 
vascular  material,  respectively  (Fig.  2).  The  boundary  r  =  Q  -  Q, 
with  outward  unit  normal  vector  n,  is  divided  into  three  mutually 
exclusive  partitions  i~y,  rq,  and  rh  for  assigning  the  Dirichlet, 
Neumann,  and  Robin  boundary  conditions,  respectively.  The  strong 
form  of  the  convection-diffusion  equation  for  the  open  domain  Q 
with  the  boundary  r  is  then  expressed  as  follows:  For  the  given 
thermal  conductivity  k  :  Q  — ►  IR2  x  IR2,  fluid  density  p  :  Qf  — >  IR, 
and  specific  heat  cp  :  Qf  — >  U  and  the  given  velocity  field 
v  :  Qf  — >  R2,  find  the  temperature  field  T  :  Q  ->  IR  such  that 

-  V  •  (kVT)  =f  in  Qs 

-  V  •  (icVT)  +  pcpY  •  VT  =f  in  Qf 

T  =  T  on  rT  (1) 

kVT  n  =  q  on  rq 

kVT  •  n  =  h(T -  T)  on  rh, 

where  /  :  Q  — >  U  denotes  the  heat  source,  T  :  />  — >  [R  the  prescribed 
temperature  along  the  Dirichlet  boundary,  and  q  :  rq  — >  U  the  ap¬ 
plied  heat  flux  along  the  Neumann  boundary,  while  h  :  rh  -►  IR 
and  Too :  rh  — >  U  describe  the  heat  transfer  coefficient  and  ambient 
temperature  along  the  convective  boundary,  respectively. 

Since  this  class  of  problems  is  characterized  by  low  Reynolds 
numbers  and  very  narrow  channels,  the  velocity  field  v  in  the 
microchannels  is  adequately  described  by  a  fully-developed  veloc¬ 
ity  profile  [25] 


Fig.  2.  Schematic  of  a  2D  domain  Q  =  Qs  u  Qf  with  the  fluid  (Qf)  and  solid  (Qs) 
phases.  The  boundary  r  with  outward  unit  normal  n  is  divided  into  regions  rT,  rq, 
and  rh  for  assigning  the  Dirichlet,  Neumann,  and  Robin  boundary  conditions, 
respectively. 
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where  rh,  A,  and  D  are  the  mass  flow  rate,  cross-sectional  area,  and 
the  diameter  of  the  microchannel,  respectively,  and  r  denotes  the 
distance  from  the  centerline  of  the  microchannel.  This  assumption 
therefore  neglects  the  transition  distance  needed  to  establish  the 
fully-developed  flow  conditions  as  the  coolant  enters  the  micro¬ 
channels,  and  the  impact  of  possible  recirculation  that  might  take 
place  in  microchannels  with  higher  curvatures  [26,27]. 

To  derive  the  weak  form  of  (1),  we  first  define  T  = 
{r  c  H1  (£2)  :  T\r?  =  7}  and  W  =  {w  c  H1  (£2)  :  w|rj.  =  0}  as  the  solu- 
tion  and  weight  function  spaces,  respectively.  The  weak  form  is 
then  written  as:  Find  TeT  such  that 

a(T,  w)  +  a(T,  w)Fh  =  ( wj )  +  (w,  q)rg  +  (w,  T00)rh  Vw  e  W,  (3) 

where  the  linear  and  bilinear  forms  in  (3)  are  given  by 


Phase  interface 


i. 

/  \  \ 

\  w 

yf  a 

j  Cf  I 

X  ^ 

/“C  j 

0* 

O  :  Standard  dofs  0  :  Generalized  dofs 


a(T,w)=  f  \7w  K\7TdQ+  f  wpcp\  VTdQ, 
J  Q  J 

a(T1w)rh=  I  hwTdr , 

Jrh 

( wj x)rh  =  (  hwTxdr , 

Jrh 

(w,q)rq=  fr  wqdr, 

(w,f)  =  [  wfdQ. 

Jq 


(4) 

(5) 

(6) 

(7) 

(8) 


The  standard  Galerkin  FEM  formulation  for  (3)  is  obtained  by  the 
discretization  of  Q  into  m  finite  elements,  Q^Qh  =  \J^1Ql  and 
employing  a  set  of  n  Lagrangian  shape  functions  {N,-(x)}"  c  Th  to 
approximate  the  solution  field  in  each  element.  Using  {N,-(x)}"  as 
the  weight  functions,  the  standard  FEM  formulation  is  then  ex¬ 
pressed  as:  Find  Th  e  Th  such  that 


a(Th,wh)  +  a(Th,w\  =  (v/,/)  -  a(whJh)  +  «q)r?  +  « 

V  wh  e  Wh,  (9) 

where 

Th(x)  =  '£Ni(x)Ti.  (10) 

Z=1 

As  discussed  in  Section  1,  an  accurate  standard  FEM  solution  for  a 
conjugate  heat  transfer  problem  requires  the  use  of  meshes  that 
conform  to  the  microchannels.  To  remove  this  cumbersome  limita¬ 
tion,  we  adopt  in  the  present  study  an  IGFEM  thermal  solver,  which, 
like  other  GFEM  formulations,  eliminates  the  need  of  using  con¬ 
forming  meshes  by  using  appropriate  enrichment  functions  to 
approximate  the  solution  field.  The  unique  feature  of  the  IGFEM, 
however,  is  how  enrichment  functions  are  evaluated  inside  the  ele¬ 
ments  cut  by  the  phase  interface  and  how  the  generalized  degrees 
of  freedom  (dofs)  are  assigned.  The  IGFEM  approximation  of  the 
solution  field  in  each  element  is  written  as  [20] 

n  nen 

T\x)  =  £N,(x)T,  +  £>,(*)«£.  (11) 

i=l  i=  1 

The  first  term  in  (1 1 )  is  the  standard  FEM  part  of  the  approximation 
with  1/  representing  the  nodal  values  of  the  temperature.  The  sec¬ 
ond  term  is  the  contribution  of  enrichment  functions  in  order  to 
capture  the  gradient  discontinuity  along  the  phase  interface,  where 
the  function  space  (i/^(x)  :  x  ->  I R}"^  is  composed  of  a  set  of  nen 
enrichment  functions,  a*  the  generalized  dofs,  and  s  is  a  scaling  fac¬ 
tor  to  avoid  sharp  gradients  in  the  enrichment  functions.  Unlike 
conventional  GFEM,  where  the  generalized  dofs  are  assigned  at 
the  nodes  of  the  original  mesh,  the  IGFEM  enrichments  are  attached 


Fig.  3.  Location  of  standard  dofs  and  generalized  dofs  at  the  interface  nodes  in  the 
IGFEM  analysis  of  a  domain  composed  of  two  phases  with  the  interface  shown  by  a 
dashed-dotted  curve. 


to  the  interface  nodes  created  at  the  intersection  points  of  elements’ 
edges  and  the  phase  interface.  Fig.  3  illustrates  the  location  of 
generalized  dofs  in  the  IGFEM  at  the  interface  nodes  in  a  domain 
composed  of  two  phases  and  discretized  with  three-node  triangular 
elements. 

A  key  advantage  of  moving  generalized  dofs  from  the  nodes  of 
original  mesh  (as  in  conventional  GFEM)  to  the  interface  nodes 
(as  in  the  IGFEM)  appears  when  assigning  Dirichlet  boundary  con¬ 
ditions  at  the  nodes  of  elements  cut  by  the  interface.  Since  the 
standard  dofs  in  the  IGFEM  directly  yield  the  nodal  values  of  the 
solution  field,  applying  the  prescribed  values  of  the  temperature 
at  these  nodes  is  similar  to  that  of  the  standard  FEM.  This  feature 
of  the  IGFEM  removes  the  need  of  techniques  such  as  the  penalty 
method  or  Lagrange  multipliers  to  assign  Dirichlet  boundary  con¬ 
ditions  at  enriched  nodes  in  conventional  GFEM  formulations 
[28,29]. 

The  enrichment  functions  in  IGFEM  are  evaluated  as  a  linear 
combination  of  Lagrangian  shape  functions  of  the  integration  ele¬ 
ments,  i.e.,  elements  created  for  the  Gauss  quadrature  of  the  parent 
element.  Fig.  4(a)  and  (b)  illustrate  two  possible  orientations  of  a 
three-node  triangular  element  with  respect  to  the  phase  interface 
and  the  corresponding  integration  elements,  which  are  given  by 

(a)  ^=Nf+Nf,  \j/2  =  Nf  +Nf\ 

(b)  ^=N«+N<2>.  1  J 

To  avoid  numerical  problems  caused  by  the  enrichment  functions  of 
elements  divided  into  integration  elements  with  high  aspect  ratios, 
we  use  the  scaling  factor  s  in  (1 1)  to  control  the  gradient  of  enrich¬ 
ment  functions.  The  problem  arises  when  an  interface  node  is  lo¬ 
cated  too  close  to  one  of  the  nodes  of  the  parent  element,  causing 
enrichment  functions  with  very  large  gradients  and  possibly  leading 
to  ill-conditioned  stiffness  matrices.  To  resolve  this  problem,  the 
enrichment  functions  are  scaled  by  the  factor  s  introduced  in  (11) 

and  defined  by  s  =  (2  nlm+i!2L2))  »  where  U  and  l2  are  the  distances 

between  the  interface  node  and  the  nodes  of  the  parent  element’s 
edge  (Fig.  5).  More  information  about  the  IGFEM,  including  the  con¬ 
vergence  study  and  applications,  can  be  found  in  [20,21]. 

3.  Validation  study 

3.1.  Experimental  set-up 

The  schematic  and  boundary  conditions  of  the  microvascular  fin 
used  for  the  validation  of  the  numerical  solver  are  illustrated  in 
Fig.  6.  The  epoxy  matrix,  with  thermal  conductivity  Ke  =  0.23  W/ 
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Fig.  4.  Evaluating  the  enrichment  functions  in  the  IGFEM.  Figs,  (a)  and  (b)  illustrate  two  scenarios  for  creating  integration  elements  and  the  associated  generalized  dofs  for 
three-node  triangular  elements  based  on  the  location  of  the  phase  interface. 


Fig.  5.  Scaling  of  the  enrichment  function  for  cases  where  the  interface  is  very  close 
to  a  node  of  the  underlying  non-conforming  mesh. 


m  K,  has  in-plane  dimensions  of  59.4  x  10.5  mm  and  a  thickness  of 
4.0  mm.  The  microvascular  fin  is  placed  on  a  60.0  x  10.0  x  0.7  mm 
copper  plate  (tcc  =  401  W/mK),  which  is  in  contact  with  a  heater 
supplied  with  constant  power  P=  2.1  W.  Other  faces  of  the  speci¬ 
men  are  in  contact  with  air  and  assumed  to  have  convective 
boundary  conditions  with  the  ambient  temperature  =  21  °C, 
and  a  heat  transfer  coefficient  h ,  to  be  determined  according  to 
the  experimental  results  obtained  in  the  absence  of  the  flow 
(Section  3.2).  The  provided  power  for  the  heater  and  hence  the  gi¬ 
ven  energy  to  the  copper  plate  is  set  such  that  the  temperature 
along  the  bottom  edge  of  the  domain  in  the  absence  of  the  coolant 
flow  is  80  °C. 

The  embedded  sinusoidal-shaped  microchannel  in  the  epoxy  fin 
with  diameter  D  =  500  pm,  amplitude  A  =  4  mm,  and  wavelength 
2=10  mm,  is  depicted  in  Fig.  6(a).  The  distance  between  the  cen¬ 
terline  of  the  microchannel  and  the  surface  over  which  the  temper¬ 
ature  field  is  recorded  by  an  infrared  camera  is  about  700  pm.  The 
coolant  flowing  in  the  microchannels  is  water  with  kw  =  0.6  W/ 
mK,  p  =  1000  kg/m3,  and  cp  =  4183  J/kg  K,  and  the  flow  rate  is 
maintained  constant  at  Q=  2  ml/min.  The  test  specimen  was  cre¬ 
ated  with  a  fugitive  organic  ink  using  the  direct  write  assembling 
[6,30]  as  shown  in  Fig.  6(c).  The  process  involves  casting  a  thin 
epoxy  substrate  upon  which  a  fugitive  ink  is  robotically  deposited 
in  a  sinusoidal  pattern  forming  a  microchannel  template.  The  pat¬ 
tern  is  then  embedded  in  epoxy  and  the  entire  sample  is  cured  at 
50  °C.  Removal  of  the  fugitive  ink  at  elevated  temperature  of 
80  °C  reveals  the  embedded  microchannel  following  the  pattern 


that  was  formed  during  deposition.  The  epoxy  used  is  Epon  862 
and  Epikure  3300  mixed  at  100:24.8  pbw.  The  fugitive  ink  used 
is  a  mixture  of  53%  microcrystalline  wax  (SP-19,  Strahl  and  Pitsch, 
Inc.)  and  47%  mineral  oil  (Fisher  Scientific).  The  ink  was  extruded 
from  a  500  pm  micronozzle  at  300  psi  at  a  printing  speed  of 
5  mm/s.  The  cured  specimen  was  cut  and  polished  to  the  desired 
dimensions  as  shown  in  Fig.  6(c).  One  side  of  the  specimen  was 
painted  matte  black  (Krylon)  for  imaging  with  an  infrared  camera. 
Propropylene  manifolds  were  attached  to  the  inlet  and  outlet 
opening  using  cyanoacrylate  glue  to  allow  access  points  for  circu¬ 
lation  of  the  coolant  through  the  microchannel. 

Active  cooling  measurements  were  accomplished  using  infrared 
imaging  of  the  microvascular  epoxy  specimen.  The  bottom  face  of 
the  epoxy  fin  was  placed  on  a  copper  plate  subjected  to  a  constant 
power  boundary  condition  resulting  in  a  constant  starting  temper¬ 
ature  in  the  absence  of  the  flow  in  the  microchannel.  A  thin  layer  of 
thermally  conductive  grease  (OmegaTherm)  was  applied  to  the 
bottom  face  of  the  composite  to  improve  thermal  contact  between 
the  copper  plate  and  the  fin.  A  resistive  heater  (Watlow,  120  V, 
15  W)  was  attached  to  the  copper  plate  and  controlled  with  a  var- 
iac  controller  (STACO  Inc.,  120  V).  A  feedback  control  thermocouple 
was  placed  on  the  top  of  the  plate  to  monitor  the  temperature  of 
the  plate  throughout  the  experiment.  A  multimeter  was  connected 
in  series  with  the  resistive  heater  to  measure  the  applied  voltage. 
The  specimen  was  allowed  to  reach  steady  state  by  waiting  for 
an  hour  after  application  of  constant  power  to  the  copper  plate. 
Fig.  6(b)  shows  a  schematic  of  the  IR  imaging  set-up  for  thermal 
management. 

Infrared  images  were  taken  using  a  DeltaTherm  1560  infrared 
camera  with  320  by  256  array  of  indium  antimonide  IR  detectors. 
Data  was  recorded  at  one  frame  per  second  using  DeltaVision 
software.  Each  data  set  contained  a  2D  temperature  field  of  the 
specimen  surface.  The  surface  temperatures  in  the  absence  of  flow 
and  steady-state  at  2  ml/min  are  depicted  in  Fig.  7(a)  and  (b), 
respectively.  The  temperature  profiles  along  lines  corresponding 
to  the  peak  (A)  and  trough  ( B )  of  the  sinusoidal  microchannel  were 
used  for  the  calibration  and  validation  of  the  2D  IGFEM  model. 
Water  was  used  as  the  coolant  and  was  introduced  into  the  sample 
at  room  temperature  (21  °C)  at  constant  flow  rate  of  2  ml/min 
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Fig.  6.  (a)  Geometry  and  boundary  conditions  of  the  validation  study  performed  on  a  thermally  loaded  microvascular  epoxy  fin  with  a  500  pm  diameter  sinusoidal  embedded 
microchannel;  (b)  schematic  of  the  experimental  set-up  with  infrared  imaging  of  the  temperature  field  on  the  surface  of  the  fin;  (c)  microvascular  epoxy  fin  specimen  with 
embedded  microvascular  channel  created  by  the  direct  ink  writing  technique. 
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Fig.  7.  IR  image  of  the  steady-state  temperature  field  on  the  surface  of  the  fin  specimen:  (a)  no  flow  case;  (b)  water  cooling  at  2  ml/min. 
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maintained  using  screw  driven  syringe  pump  (KD  scientific,  Model 
210). 

3.2.  Calibration  and  validation  of  the  IGFEM  thermal  solver 

The  calibration  component  of  this  study  consists  in  evaluating 
the  heat  transfer  coefficient  h ,  associated  with  the  surfaces  of  the 
specimen  in  contact  with  air,  which  is  needed  to  determine  the 
dissipated  heat  from  these  surfaces.  To  that  effect,  we  use  the  data 
obtained  with  IR  measurements  in  the  absence  of  the  flow,  as 
shown  in  Fig.  7(a).  Since  a  2D  model  is  adopted  for  approximating 
the  temperature  field,  we  incorporate  the  dissipated  heat  from  the 
faces  of  the  specimen  perpendicular  to  the  computational  domain 
by  adopting  the  classical  fin  equation  expressed  as 

-V-(kVT)  =f-f(T-Tx),  (13) 

where  P  and  A  are  the  perimeter  and  area  of  the  face  of  the  speci¬ 
men  in  contact  with  the  copper  plate,  respectively. 

In  the  absence  of  the  flow,  the  temperature  field  is  quasi  ID  as 
apparent  both  experimentally  (Fig.  7(a))  and  numerically 
(Fig.  8(a)).  The  temperature  variation  in  the  fin  along  a  line  perpen¬ 
dicular  to  the  copper  plate  can  be  analytically  expressed  as 

T  -  Tx  _  cosh  P{L  - y)  +  ^  sinh  fl{L  - y) 

T -Too  cosh +  ±  sinh [SL 

where  p  =  with  hf  =  ^  denoting  the  fin  coefficient.  The  bound¬ 
ary  conditions  consist  in  a  prescribed  temperature  T  =  80  °C  along 
the  bottom  edge  of  the  domain  and  convective  boundary  conditions 
along  the  top  edge.  A  direct  comparison  between  this  relation  and 
the  experimental  results  obtained  along  line  segment  A  in 


Fig.  7(a)  yields  a  least  square  fit  value  of  the  film  coefficient  h  equal 
to  4.9  W/m  K.  A  comparison  between  analytical,  numerical,  and 
experimental  results  is  presented  in  Fig.  8(b),  showing  excellent 
agreement. 

The  next  step  is  to  evaluate  the  equivalent  flow  rate  for  the  2D 
solver  that  yields  a  similar  temperature  field  as  that  obtained 
experimentally.  The  experimental  results  for  the  temperature  field 
after  steady  state  are  shown  in  Fig.  7(b).  In  this  case,  the  copper 
plate  average  temperature  is  reduced  to  65  °C  and  due  to  the  high 
thermal  conductivity  of  the  copper,  this  temperature  is  almost  con¬ 
stant  at  every  point  along  the  length  of  the  plate.  The  equivalent 
flow  rate  for  the  2D  model  is  then  evaluated  such  that  it  yields 
the  least  square  fit  of  the  temperature  profiles  from  the  experimen¬ 
tal  and  numerical  results  along  Sections  A  and  B,  shown  in  Fig.  7(b). 
This  calculation  gives  Qeq  =  0.724  ml/min,  for  which  the  steady  state 
temperature  field  obtained  from  the  IGFEM  solution  and  the  tem¬ 
perature  profiles  along  Sections  A  and  B  are  presented  in  Fig.  9. 

To  validate  the  2D  IGFEM  solver,  we  need  to  examine  the  per¬ 
formance  of  this  method  for  evaluating  the  temperature  field  in  a 
microvascular  material  with  similar  dimensions  and  boundary 
conditions,  but  with  a  different  configuration  of  the  embedded 
microchannel  as  illustrated  in  Fig.  10(a).  The  heat  transfer  and  fin 
coefficients  for  this  specimen,  which  are  functions  of  the  dimen¬ 
sions  of  the  domain,  are  similar  to  those  of  the  specimen  employed 
in  the  calibration  study  (Fig.  6(c)).  For  the  embedded  microchannel 
shown  in  Fig.  10(a),  the  wavelength  of  the  channel  is  reduced  to 
2  =  5  mm,  and  a  flow  rate  of  Q.  =  2  ml/min  is  used  for  active  cooling 
of  the  material.  The  recorded  temperature  field  with  the  IR  camera 
over  the  surface  of  the  actively-cooled  specimen,  adjacent  to  the 
embedded  microchannel,  is  depicted  in  Fig.  10(b). 

The  temperature  field  evaluated  by  the  2D  IGFEM  solver  for  the 
validation  problem  is  illustrated  in  Fig.  11(a),  where  the  equivalent 
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Fig.  8.  (a)  Temperature  field  obtained  with  the  2D  IGFEM  solver  in  the  absence  of  the  flow;  (b)  experimental,  numerical  and  analytical  temperature  profiles  along  Section  A.  A 
heat  transfer  coefficient  of  h  =  4.9  W/mK  is  employed  to  evaluate  the  analytical  and  numerical  solutions. 
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Fig.  9.  Calibration  of  the  IGFEM  solver  using  the  thermally  loaded  domain  with  an  embedded  microchannel  with  2  =  10  mm:  (a)  steady-state  temperature  field  of  the 
actively-cooled  fin  at  2  ml/min  flow  rate  obtained  with  the  IGFEM;  (b,c)  experimental  and  numerical  temperature  profiles  along  Sections  A  and  B,  respectively. 
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Fig.  10.  (a)  Microvascular  fin  specimen  with  an  embedded  microchannel  with  wavelength  2  =  5  mm  and  amplitude  A  =  4  mm;  (b)  IR  image  of  the  steady-state  temperature 
field  after  active  cooling  at  2  ml/min  flow  rate. 


flow  rate  obtained  from  the  calibration  study  ( Qeq  =  0.724  ml/min) 
is  implemented  in  the  numerical  solution.  The  temperature 
profiles  along  Sections  A  and  B  (shown  in  Fig.  10(b))  are  presented 
in  Fig.  11(b)  and  (c),  respectively,  and  compared  with  experimental 
values  at  similar  locations.  Accoring  to  this  figure,  the  2D  IGFEM 
solution  has  a  good  agreement  with  experimental  data,  especially 
in  predicting  the  temperature  along  the  top  edge  of  the  domain. 
Considering  all  the  simplifying  assumptions  implemented  in  the 
2D  computational  model,  this  study  proves  the  reliability  of  the 
IGFEM  in  approximating  the  temperature  field  in  the  actively- 
cooled  microvascular  fin. 


4.  Shape  optimization  of  the  sinusoidal  microchannel 

In  this  section,  we  aim  at  determining  the  optimal  configuration 
of  the  embedded  sinusoidal-shaped  microchannel  in  the  epoxy  fin 
shown  in  Fig.  12.  The  height  and  thickness  of  this  domain  are 
similar  to  those  of  the  microvascular  fin  specimen  depicted  in 
Fig.  6  used  for  the  calibration  and  validation  of  the  IGFEM  thermal 
solver.  The  heat  transfer  and  fin  coefficients  for  the  shape  optimiza¬ 
tion  problem  are  similar  to  those  obtained  in  the  previous  section. 
The  length  of  the  microvascular  fin,  however,  has  been  increased 
to  180  mm  to  investigate  the  effect  of  the  specimen  length  on  the 
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Fig.  11.  Validation  of  the  IGFEM  solver  based  on  the  microvascular  fin  with  -1  =  5  mm:  (a)  steady-state  temperature  field  of  the  actively-cooled  fin  at  2  ml/min  flow  rate 
obtained  with  the  IGFEM;  (b,c)  experimental  and  numerical  temperature  profiles  along  Sections  A  and  B,  respectively. 
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Fig.  12.  Schematic  and  boundary  conditions  of  the  microvascular  fin  investigated  in  the  shape  optimization  study. 


optimal  shape  of  the  embedded  microchannel.  The  thermal  bound¬ 
ary  conditions  consist  of  a  convective  boundary  condition  along  the 
left  and  right  sides  of  the  2D  computational  domain  and  constant 
heat  flux  q  =  2000  W/m2  along  the  bottom  edge.  Two  different 
conditions  are  considered  along  the  top  edge  to  study  their  impact 
on  the  predicted  optimal  configuration  of  the  embedded  micro- 
channel:  (i)  a  fixed  temperature  of  T  =  20  °C  and  (ii)  a  convective 
boundary  condition  with  h  =  4.9  W/mK  and  T  =  21  °C. 

The  optimization  problem  aims  at  determining  the  optimal 
amplitude,  A ,  and  wavelength,  2,  of  the  sinusoidal  microchannel 
that  minimizes  a  set  of  three  objective  functions  <D  =  {0,-}f=1:  (i) 
the  maximum  temperature  of  the  polymer  matrix,  (ii)  its  void  vol¬ 
ume  fraction,  and  (iii)  the  pressure  drop.  Details  on  the  evaluation 
of  these  three  objective  functions  are  provided  below.  The  entrance 
of  the  microchannel  is  at  the  left  edge  and  the  minimum  distance 
between  its  centerline  and  the  bottom  edge  of  the  domain  is 
580  pm  for  all  configurations.  Similar  to  the  previous  section, 
water  is  employed  as  the  coolant  with  Q.  =  1  ml/min,  which,  based 
on  the  study  presented  in  the  previous  section,  corresponds  to  an 
equivalent  flow  rate  of  Qeq  =  0.362  ml/min  for  the  2D  model. 

The  minimum  and  maximum  values  of  the  wavelength  2  are  set 
to  7.5  and  30  mm,  respectively,  with  increments  of  2.5  mm  while 
the  minimum  and  maximum  values  of  A  are  0  and  4.8  mm,  respec¬ 
tively,  with  increments  of  400  pm.  The  aforementioned  parameters 
provides  a  total  number  of  121  different  configurations  for  which 
the  small  size  of  the  search  space  allows  for  analyzing  and 
evaluating  the  objective  functions  for  all  possible  cases.  This  is 


specifically  possible  by  using  the  IGFEM  solver,  for  which  a  single 
non-conforming  structured  mesh  built  on  a  1080  x  63  grid  is  em¬ 
ployed  to  evaluate  the  temperature  field  associated  with  all  the 
microchannel  configurations.  The  objective  functions  considered 
in  this  optimization  problem  are  as  follows: 

Maximum  temperature :  Reducing  the  maximum  temperature  of 
the  domain  is  the  main  goal  of  designing  actively-cooled  microvas¬ 
cular  systems.  We  thus  define  the  first  objective  function  as 

(f>i  =  Tmax  :  Qh  -►  |R  :=  ma xTh,  (15) 

where  t1  is  obtained  from  (9)  using  the  IGFEM  solver.  For  this  par¬ 
ticular  optimization  problem,  the  maximum  temperature  is  typi¬ 
cally  obtained  along  the  lower  edge  of  the  specimen,  along  which 
the  heat  flux  is  applied. 

Void  volume  fraction :  While  reducing  the  maximum  tempera¬ 
ture  of  the  domain  constitutes  the  primary  objective,  a  successful 
design  for  the  microvascular  network  must  be  minimally  invasive 
with  respect  to  the  mechanical  properties  of  the  material,  and  in 
particular  its  strength  and  stiffness.  As  a  first  approximation,  the 
impact  of  the  microchannel  on  these  mechanical  properties  can 
be  related  to  the  void  volume  fraction  associated  with  the  embed¬ 
ded  network.  The  second  objective  function  thus  consists  in 
minimizing 

fj  nD2ds 

(t>2  =  Vf  <eR  :=  - ,  (16) 
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where  VQh  is  the  total  volume  of  the  microvascular  material,  Lf  and 
D  are  respectively  the  length  and  diameter  of  the  channel,  and  ds  is 
a  longitudinal  element  at  the  centerline  of  the  microchannel.  It 
must  be  noted  that,  due  to  the  constant  diameter  of  the  microchan¬ 
nels,  minimizing  this  objective  function  corresponds  to  minimizing 
the  length  of  the  embedded  microchannel  for  the  current  system. 

Pressure  drop :  The  third  objective  function  is  defined  as  the 
required  pressure  to  move  the  coolant  between  the  inlet  and  the 
outlet  of  the  microchannel: 

03  =  Ap  6  R  :=  max  | pin  - p0,„|.  (17) 

This  pressure  drop  is  directly  related  to  the  flow  efficiency  of  the 
microchannel,  because  the  power  required  to  circulate  the  fluid  is 
expressed  as  P  =  mAp/p.  Assuming  a  fully-developed  flow  velocity 
profile  as  described  by  (2),  the  pressure  drop  in  the  microchannel  is 
obtained  from  the  Hagen-Poiseuille  law  [25] 

(18) 

JLf  pnD 

where  p  is  the  dynamic  viscosity  of  the  coolant,  which  is  a  decreas¬ 
ing  function  of  its  temperature  Tf.  In  the  current  study,  we  use  the 
Seeton  fit  [31]  to  approximate  the  temperature  dependence  of  the 
viscosity  for  water: 

p{Tf)  =  Cx  lO7^,  (19) 

where  C=  2.414  x  10-5  N  s/m2,  7)  =  247.8  K,  and  T2  =  140  K.  The  dy- 
namic  viscosity  of  water  is  then  evaluated  along  each  longitudinal 
cross  section  ds  of  the  microchannel  according  to  (19)  to  determine 
the  pressure  drop  in  that  cross  section. 


4.1.  Prescribed  temperature  along  the  top  edge 

Assuming  a  prescribed  temperature  of  Tt  =  20  °C  along  the  upper 
edge  of  the  domain  (Fig.  12),  the  numerical  predictions  obtained 
with  the  IGFEM  thermal  solver  are  presented  in  Fig.  13.  The  values 
of  the  objective  functions  for  each  individual,  i.e.,  each  data  point 
corresponding  to  one  of  the  configurations  of  the  microchannel, 
are  depicted  in  Fig.  13(a).  The  void  volume  fraction  versus  the  max¬ 
imum  temperature  of  the  domain  for  each  individual  is  shown  in 
Fig.  13(b)  to  better  illustrate  the  variations  of  these  two  objective 
functions  for  different  configurations.  As  shown  there,  the  individ¬ 
uals  corresponding  to  the  lowest  values  of  the  void  volume  fraction 
maintain  a  relatively  low  maximum  temperature  in  the  domain. 
The  maximum  temperature  of  the  domain  then  increases  with 
the  void  volume  fraction  until  it  reaches  a  peak  value  beyond 
which,  the  microchannels  with  higher  void  volume  fractions  are 
more  effective  in  reducing  the  maximum  temperature. 

The  variations  of  the  maximum  temperature  of  the  material 
with  respect  to  the  wavelength  and  the  amplitude  of  the  micro¬ 
channels,  and  with  respect  to  its  amplitude  for  X  =  7.5  mm  are 
shown  in  Fig.  13(c)  and  (d),  respectively.  As  shown  in  Fig.  13(d),  a 
straight  microchannel  between  the  inlet  and  the  outlet  (A  =  0 
mm)  reduces  the  maximum  temperature  of  the  domain  to 
Tmax  =  60.8  °C.  Increasing  the  amplitude  of  the  microchannel  up  to 
A  =  3.6  mm  not  only  increases  the  void  volume  fraction,  but  also 
is  less  effective  in  decreasing  the  maximum  temperature  of  the 
domain.  After  this  point,  however,  microchannels  with  higher 
amplitudes  yield  more  effective  configurations  for  reducing  the 
maximum  temperature  compared  to  the  straight  channel,  although 
this  comes  at  the  price  of  a  higher  void  volume  fraction  and 
pressure  drop. 


Fig.  13.  Optimization  results  for  the  problem  shown  in  Fig.  12  with  a  prescribed  temperature  along  the  top  edge:  (a)  representation  of  three  and  (b)  two  objective  functions 
for  each  individual;  (c)  maximum  temperature  versus  the  wavelength  k  and  microchannel  amplitude  A\  (d)  maximum  temperature  in  the  fin  versus  the  microchannel 
amplitude  for  l  =  7.5  mm. 
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Fig.  14.  Temperature  field  associated  with  the  configuration  labeled  as  (a-c)  in  Fig.  13(b). 
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Fig.  15.  Temperature  profile  along  the  bottom  edge  of  domain  associated  with  the  microchannels  shown  in  Fig.  14. 


The  temperature  fields  for  the  individuals  marked  as  (a-c)  in 
Fig.  13(b)  are  illustrated  in  Fig.  14(a-c),  respectively.  As  expected, 
the  temperature  field  associated  with  the  straight  microchannel 
shown  in  Fig.  14(a)  corresponds  to  the  best  configuration  for  min¬ 
imizing  the  void  volume  fraction  and  pressure  drop.  The  micro- 
channel  shown  in  Fig.  14(b)  with  the  minimum  wavelength  and 
maximum  amplitude  is  the  most  effective  for  reducing  the  maxi¬ 
mum  temperature,  although  this  configuration  yields  the  highest 
void  volume  fraction  and  pressure  drop.  Finally,  the  microchannel 
corresponding  to  the  highest  temperature  over  the  domain  and  the 
associated  temperature  field  are  depicted  in  Fig.  14(c). 

The  impact  of  the  microchannel  configuration  on  reducing  the 
maximum  temperature  of  the  domain  can  be  better  understood 
by  studying  Fig.  15.  This  figure  illustrates  the  temperature  profile 
along  the  bottom  edge  of  the  domain,  where  the  maximum  tem¬ 
perature  occurs,  for  the  three  microchannels  shown  in  Fig.  14. 
The  temperature  profile  associated  with  the  straight  microchannel 
has  a  roughly  constant  rate  of  increase  along  this  edge.  For  this 
microchannel,  the  only  mechanism  that  contributes  to  active  cool¬ 
ing  of  the  domain  is  removing  the  heat  through  the  outgoing  flow. 
Therefore,  as  we  move  towards  the  right  edge  of  the  domain,  the 
maximum  temperature  of  the  material  increases  due  to  the  in¬ 
crease  of  the  coolant  temperature  and  hence  the  reduction  of  its 
cooling  capacity. 

The  temperature  profile  along  the  bottom  edge  of  the  domain 
associated  with  the  sinusoidal  channel  depicted  in  Fig.  14(b) 
(A  =  4.8  mm,  X  =  7.5  mm)  reveals  a  different  thermal  response.  As 
shown  in  Fig.  15,  after  only  a  few  oscillations  of  the  microchannel, 


the  temperature  along  the  bottom  edge  reaches  a  periodic  state 
where  the  inflow  and  outflow  temperatures  of  the  fluid  at  the 
beginning  and  the  end  of  each  phase  are  similar.  This  is  achieved 
through  the  heat  exchange  between  the  lower  and  upper  regions 
of  the  domain,  i.e.,  the  warmer  and  colder  areas,  respectively,  such 
that  the  amount  of  the  extracted  and  transmitted  heat  in  each 
region  are  equal.  As  explained  in  Section  1,  configurations  that 
can  effectively  redistribute  the  heat  inside  the  domain  are  better 
options  for  active  cooling  of  domains  similar  to  that  of  the  current 
problem  with  large  distances  between  the  inflow  and  outflow 
points.  Compared  to  the  straight  channel,  Fig.  15  verifies  that  the 
sinusoidal  microchannel  with  A  =  4.8  and  A  =  7.5  mm  is  more  effec¬ 
tive  in  reducing  the  maximum  temperature  of  the  material  if  the 
length  of  the  domain  exceeds  10  cm.  Moreover,  the  sinusoidal  con¬ 
figuration  maintains  the  same  maximum  temperature  regardless 
of  the  length  of  the  domain  while  the  maximum  temperature  asso¬ 
ciated  with  the  straight  microchannel  has  an  increasing  rate  with 
increasing  the  length  of  the  domain. 

The  temperature  profile  associated  with  the  microchannel  cor¬ 
responding  to  the  highest  maximum  temperature  over  the  domain 
with  A  =  3.6  and  A  =  30  mm  shows  a  similar  behavior  as  the  previ¬ 
ous  case.  However,  the  large  wavelength  and  average  amplitude  of 
this  channel  which  do  not  let  effective  heat  exchange  with  the  cold 
region  near  the  upper  edge  of  the  domain  yields  a  much  higher 
maximum  temperature  and  has  not  reached  the  periodic  state 
within  the  length  of  the  domain.  It  must  be  noted  that  as  shown 
in  Fig.  15,  the  increasing  rate  of  the  temperature  along  the  bottom 
edge  associated  with  this  sinusoidal  configuration  is  lower  than 
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that  of  the  straight  microchannel.  Thus,  it  is  not  difficult  to  predict 
that  as  the  length  of  the  domain  increases,  even  this  configuration 
will  eventually  be  more  effective  for  active  cooling  of  the  domain 
compared  to  the  straight  channel. 


42.  Convective  boundary  conditions  along  the  top  edge 

Fig.  16  illustrates  the  optimization  results  corresponding  to 
those  shown  in  Fig.  13,  but  assuming  convective  boundary  condi¬ 
tions  along  the  top  edge  of  the  domain.  As  shown  there,  the  change 
of  boundary  conditions  has  a  substantial  effect  on  the  optimal 
configuration  of  the  microchannel,  as  the  straight  channel  repre¬ 
sents  the  optimal  shape  with  respect  to  all  the  objective  functions 
in  this  case.  In  other  words,  increasing  either  the  wavelength  or  the 


amplitude  of  the  microchannels  not  only  exacerbates  the  void 
volume  fraction  and  the  pressure  drop,  but  also  has  an  adverse  im¬ 
pact  on  the  maximum  temperature  of  the  material. 

The  temperature  fields  and  microchannel  configurations  for  the 
three  individuals  labeled  as  (a-c)  in  Fig.  16(b)  are  depicted  in 
Fig.  17(a-c),  respectively.  The  corresponding  temperature  profiles 
along  the  bottom  edge  of  the  fin  are  illustrated  in  Fig.  18.  The  tem¬ 
perature  field  associated  with  the  straight  channel,  depicted  in 
Fig.  17(a),  is  very  similar  to  that  shown  in  Fig.  14(a)  with 
prescribed  temperature  along  the  top  edge.  Comparing  Figs.  18 
and  15,  we  observe  that  the  variations  of  the  temperature  along 
the  bottom  edge  for  both  cases  are  almost  identical.  However,  un¬ 
like  the  problem  with  prescribed  temperature  along  the  top  edge, 
the  straight  microchannel  is  the  best  configuration  for  reducing 
the  maximum  temperature  of  the  domain  in  the  current  case. 


vf 


(b) 


Fig.  16.  Optimization  results  for  the  problem  shown  in  Fig.  12  with  convective  boundary  conditions  along  the  top  edge:  (a)  representation  of  three  and  (b)  two  objective 
functions  for  each  individual;  (c)  maximum  temperature  versus  the  wavelength  X  and  microchannel  amplitude  A\  (d)  maximum  temperature  in  the  fin  versus  the 
microchannel  amplitude  for  X  =  7.5  mm. 
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Fig.  17.  Temperature  field  associated  with  the  individuals  labeled  as  (a-c)  in  Fig.  17(a). 
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Fig.  18.  Temperature  profile  along  the  bottom  edge  of  domain  associated  with  the  microchannels  shown  in  Fig.  17. 


Fig.  17(b)  and  (c)  illustrate  the  microchannels  associated  with 
the  highest  void  volume  fraction  (and  pressure  drop)  and  highest 
maximum  temperature  over  the  domain,  respectively.  As  shown 
in  Fig.  18,  despite  the  sinusoidal  shape  of  the  microchannel  and 
thus  the  heat  exchange  between  the  lower  and  upper  edges  of 
the  domain,  the  temperature  along  the  bottom  edge  never  shows 
a  periodic  behavior.  The  reason  is  that  when  the  coolant  flows 
towards  the  top  edge,  it  increases  the  temperature  of  the  matrix 
in  that  region  and  hence  the  oscillations  of  the  microchannel  for 
exchanging  the  heat  between  the  top  and  bottom  edges  become 
ineffective  for  reducing  the  temperature  of  the  coolant  along  the 
length  of  the  domain.  This  behavior  is  completely  different  from 
the  previous  case  with  the  prescribed  temperature  along  the  top 
edge,  where  the  rate  of  increase  of  the  temperature  was  descend¬ 
ing  and  a  periodic  state  was  eventually  emerging.  Thus,  the  only 
efficient  configuration  for  active  cooling  of  the  problem  shown  in 
Fig.  12  with  convective  boundary  conditions  along  the  top  edge 
is  a  straight  microchannel  between  the  inlet  and  outlet.  However, 
the  straight  configuration  is  only  effective  along  a  short  length  due 
to  the  increasing  temperature  of  the  coolant,  which  deteriorates  its 
cooling  capacity.  For  longer  domains  with  similar  boundary  condi¬ 
tions,  multiple  inlets  and  outlets  must  be  considered  at  appropriate 
distances  along  the  length  of  the  material  to  avoid  excessive 
increase  of  the  temperature  of  the  coolant  before  existing  the 
domain. 

5.  Conclusions 

The  optimal  configuration  of  the  embedded  sinusoidal  micro- 
channel  in  an  epoxy  fin  was  studied.  An  interface-enriched  GFEM 
(IGFEM)  was  employed  to  evaluate  the  temperature  field  for  differ¬ 
ent  configurations  of  the  microchannel.  This  method  provides 
independence  between  the  structure  of  the  finite  element  mesh 
and  geometry  of  the  problem.  Moreover,  the  IGFEM  provides 
several  advantages  over  conventional  GFEM  including  the  easy 
assignment  of  Dirichlet  boundary  conditions  at  enriched  elements 
and  low  computational  cost.  The  IGFEM  solver  was  calibrated  and 
validated  using  the  experimental  data  obtained  from  thermal  tests 
performed  on  a  microvascular  epoxy  fin.  The  heat  transfer  and  fin 
coefficients  of  the  domain  were  evaluated  using  the  results  of  the 
test  in  the  absence  of  active  cooling,  while  the  data  obtained  from 
the  actively-cooled  fin  were  employed  to  extract  the  equivalent  2D 
flow  rate.  The  IGFEM  model  was  then  validated  by  evaluating  the 
temperature  field  in  a  similar  microvascular  fin  with  a  different 
microchannel  configuration  and  comparing  that  with  the  experi¬ 
ment.  Finally,  we  employed  the  IGFEM  solver  to  determine  the 
optimal  microchannel  configurations  by  finding  the  wavelength 
and  amplitude  of  the  sinusoidal  channel  that  minimize  the 


maximum  temperature,  void  volume  fraction,  and  the  pressure 
drop  for  circulating  the  fluid.  We  also  studied  the  impact  of  bound¬ 
ary  conditions  on  the  optimal  configuration  of  the  channels,  which 
proved  to  be  a  decisive  factor  in  the  design  of  the  microvascular 
fin. 
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